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Abstract 



We analyze in full mathematical rigor the most general quartically perturbed in- 
i-pn ', variant probability measure for a random tensor. Using a version of the Loop Vertex 

Expansion (which we call the mixed expansion) we show that the cumulants write 
' ■^ . as explicit series in 1/N plus bounded rest terms. The mixed expansion recasts the 

problem of determining the subleading corrections in 1/N into a simple combinatorial 
problem of counting trees decorated by a finite number of loop edges. 

As an aside, we use the mixed expansion to show that the (divergent) perturbative 
expansion of the tensor models is Borel summable and to prove that the cumulants 
respect an uniform scaling bound. In particular the quartically perturbed measures 
fall, in the N — > oo limit, in the universality class of Gaussian tensor models. 

vo' 

tj- ! 1 Introduction 

o 

Tensor models [T] generalize matrix models [2j [3] to higher dimensions and provide 
the analytical tool for the study of random geometries in three and more dimensions. 
Matrix models are probability distributions for random matrices and their moments 
can be evaluated via a perturbative expansion in ribbon Feynman graphs representing 
surfaces [3]. A crucial tool in matrix models is the 1/N expansion discovered by 't 
Hooft [3]. The perturbative series of matrix models can be reorganized as a series 
in 1/N (where N is the size of the matrix) indexed by the genus. At leading order 
only planar graphs [5] contribute and a matrix model undergoes a phase transition to 
a theory of random continuum surfaces when the coupling constant is tuned to some 
critical value [7] . 

However many questions concerning the 1/N expansion of matrix models remain 
unanswered. First and foremost the 1/N expansion seems inherently a perturbative 
tool: in order to perform it one first performs the expansion in the coupling. The 
perturbative series is not summable hence it is not clear if the subsequent steps are 
mathematically meaningful. Furthermore, matrix models are defined for a fixed sign of 
the coupling constant (say positive) for which the perturbation is stable. However the 
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phase transition to continuum surfaces takes place when tuning the coupling constant to 
a negative critical value. This tuning is meaningful after restricting to the leading order 
planar series (which is absolutely convergent for both signs of the coupling constant), 
but what, if any, is the meaning of this tuning to criticality beyond perturbation theory? 
Can one actually reach this phase transition by an analytic continuation starting from 
a matrix model? These questions have not yet been satisfactory answered. 

Matrix models have been generalized in higher dimensions to tensor models and 
group field theories P [9l [lOl HH H21 HH HH [15]. They encode models of random 
geometries relevant for quantum gravity. However, progress has been slow (beyond 
model building) in higher dimensions mainly due to the lack of a 1/N expansion for 
tensor models. 

This has changed with the advent of the 1/N expansion, initially for the colored 
[161 fT7] models, and subsequently for all invariant tensor models [18J . Indeed the 
perturbative series of tensor models supports a 1/N expansion [191 EOj [211 1221 [23] 
indexed by the degree, a positive integer which plays in higher dimensions the role of 
the genuqj. The leading order melonic [24[ [25] graphs triangulate the D-dimensional 
sphere in any dimension |19 [ I20 [ |2T] and, like their two dimensional counterparts, tensor 
models undergo a phase transition to a theory of continuous random spaces when tuning 
to criticality. These results led to significant progress in the understanding of random 
geometries in higher dimensions [261 123 [281 (23 EQl ED E21 (Ml ESI EH ED [38] and the 
related critical phenomena. Tensor models have been generalized to renormalizable 
(and generically asymptotically free) tensor field theories, [39 ] HO ] HT ] B2 ] [33 ] H^ | H5 ] H6] 
leading to the formulation of the "tensor track" approach to quantum gravity |47[ [48] . 

Tensor models have been shown to exhibit a powerful universality property |49j : 
all invariant models respecting a uniform bound on the cumulants become Gaussian in 
the N — > oo limit, but the covariance of the limiting Gaussian is a non trivial function 
of the parameters of the model. 

The same questions concerning the status of the 1/N for matrix models can also be 
formulated for tensor models. Can one give a meaning to the 1/N expansion beyond 
perturbation theory? Can one reach the phase transition point by analytic continua- 
tion? This paper answers the first question and establishes the 1/N expansion in the 
constructive field theory sense |50] . In order to achieve this result we consider the most 
general tensor model with a relevant quartic interactional. We generalize to tensors the 
Loop Vertex Expansion (LVE) [51[ [52] initially introduced for matrices_) and analyze it 
in detail. The LVE expresses any cumulant as an absolutely convergent series indexed 
by trees. In this series the coupling constant A appears in two places: first as an overall 
factor for each tree and second in the contribution of each tree. This allows for instance 
to prove that, like in matrix models, the perturbative series of tensor models is Borel 
summable in A uniformly in N. 

However, in the case of tensors the LVE is much more powerful. Although rather 
involved at first sight, this formulation has two important advantages. First, in the 
contribution of each tree, the coupling constant always appears rescaled as \/N D ~ l 
where D is the rank of the tensor. Performing a Taylor expansion in X/N one 

proves that the rest term is suppressed in powers of \\\/N D ~ 2 . We call this expansion 



Unlike the genus, the degree is not a topological invariant. 

2 Othcr quartic interactions can be added but are suppressed in 1/N. 
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the mixed expansion. The mixed expansion provides the 1/N expansion of the model 
beyond perturbation theory: although shifted from iV D_1 to N D ~ 2 , the scaling with N 
suppresses the rest term as long as D > 3. This does not hold for the case of matrices, 
D = 2. Second, at each order in 1/N, one must sum the series indexed by trees and this 
can be done explicitly. Somewhat surprisingly, the mixed expansion turns out to be 
the appropriate computational tool for the study of the 1/N series. It reorganizes the 
corrections in 1/N in terms of trees with a finite number of loop edges, and counting 
such trees is a straightforward combinatorial problem. 

The critical behavior of each order in 1/N is governed order by order by the diver- 
gence of the series indexed by trees. In particular this proves that all terms will diverge 
for the same critical constant. The mixed expansion (more precisely the 1/N expansion 
derived from it) is an explicitly perturbation in 1/N around the leading order theory 
of trees with no loop edges. The leading order (melonic) theory acts effectively as a 
new "vacuum" around which the subleading terms in 1/N act as perturbations. 

Beyond providing a tool to analyze order by order the subleading behavior in 1/N 
of the tensor models and an avenue towards establishing their double scaling limits, 
the present study is a needed step in order to analytically continue the cumulants to 
the critical constant of the phase transition to the continuum theory. 

The situation is more subtle for matrix models. The rest terms in the mixed ex- 
pansion are not suppressed in powers of 1/N. In order to obtain a non perturbative 
definition of the 1/N series for matrices one needs to refine further the mixed expansion 
introduced in the present paper. We believe that an in depth study of the corrections in 
1/N in the case of tensors will provide the guide towards obtaining the non perturbative 
1/N expansion for matrices. 

We emphasize that the techniques of this paper rely heavily on the LVE formalism: 
some familiarity with |51} I52| would greatly benefit the reader. 

This paper is organized as follows. Section [2] is an introductory section. In it we 
present in some details the framework of tensor models, we briefly recall the notion of 
Borel summability and we introduce at length the various notions used in the sequel. 
In section [3] we introduce the most general model with a relevant quartic perturbation, 
and we state and comment our main theorems. Section [J] contains the proofs of our 
results. 

2 Prerequisites 

In this section we present the general setting of random tensor models and introduce 
the various notions and notations we will use in the sequel. 

2.1 Generalities on Tensor Models 

We start by a brief overview of the general framework of invariant tensor models pre- 
sented in detail in [18|, I49j . This allows us to introduce some relevant notions and 
notations and to state the universality theorem for random tensors. Later on we will 
prove that the quartically perturbed measure we deal with in this paper obeys indeed 
the universality theorem. 

We consider rank D covariant tensors T n i n o, with n 1 ,n 2 ,...n £ {1,...N}. 
having no symmetry under permutation of their indices. The tensors transform 



under the external tensor product of D fundamental representations of the unitary 
group U(N) (that is the unitary group acts independently on each index). The complex 
conjugate tensor, T n i iiin D is a rank D contravariant tensor. A tensor and its complex 
conjugate can be seen as collections of N complex numbers supplemented by the 
requirement of covariance under base change. 



T„i 



2_j U a 1 n 1 ■ ■ ■ V a D n D ^n 1 ...n D 



^a 1 ...a D - 2^ U n i n i ...V E D n DT n i M D 



(1) 



where the indices of the complex conjugated tensor are denoted conventionally with a 
bar. We emphasize that the unitary operators U, . . .V are all independent We denote 
n the D-tuple of integers (n , . . . n ), and we take D > 3. 

Any invariant polynomial in the tensor entries can be expressed in terms of the 
trace invariants built by contracting in all possible ways pairs of covariant and con- 
travariant indices in a product of tensor entries (see |54| [55] for a direct proof relying 
on averaging over the unitary group). The trace invariants are one to one with closed 
D-colored graphs. 

Definition 1. A bipartite closed L>-colored graph is a graph B = (V(B),£(B)) 
with vertex set V(B) and edge set £{B) such that: 

• V(B) is bipartite, i.e. the vertex set writes as V(B) = A(B) U A(B), such that 
V7 G £(B), then I = (v, v) with v G A(B) and v S A(B). Their cardinalities satisfy 
\V(B)\ = 2\A(B)\ = 2\A(B)\. We call v e A(B) the white vertices and v G A{B) 
the black vertices of B. 

• The edge set is partitioned into D subsets £(B) = [J i= i£ l (B), where £ 1 {B) = 
{P = (v, v)} is the subset of edges with color i. 



• All vertices are D-valent with all edges incident to a given vertex having distinct 
colors. 

The graph associated to an invariant (see Figure [1] for some examples) is obtained 
as follows. We represent every T n i n o (respectively T n i_ n r>) by a white vertex v 
(respectively a black vertex v). The position of an index becomes a color, n 1 has color 
1, n 2 has color 2 and so on. We represent by an edge the contraction of an index ri l on 
T n i n D with an index n 1 of T fi i___ fl D. The edges l l = (v,v) G £ 1 {B) inherit the color i 
of the indices and always connect a white and a black vertex. 




T, 
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Figure 1: .D-colored graphs representing trace invariants. 



One can consider more generally tensors transforming under the external tensor product of fundamental 
representations of unitary groups of different sizes U{N 1 ) M U(N 2 ) E • • • El U(N D ) with N t ^ Nj. 



The trace invariant associated to B is 

D 

Tr B (T,t)=^c n T^t, o , c=n n **<*<■ (2) 

n,n v,veV{B) «=1 l l = (v,v)e£ i (B) 

The trace invariant associated to B factors over its connected components. We 
call the invariant connected if B is connected. We denote k(B) the number of white 
(or black) vertices of the graph B and C(B) the number of connected components of 
B, which we label B p . Colored graphs are dual to D-dimensional abstract simplicial 
pseudo-manifolds [161 El EH [57] . 

A random tensor is a collection of N complex random variables whose joint distri- 
bution is encoded in the cumulants (connected moments) of k entries T and k entries T, 
K 2k [Tfj, , T=>_ . . . T=* , T* 1. We consider only even distributions, that is the cumulants 

1 ™ k _ 

are nontrivial only if k = k. 

Definition 2. The probability distribution [In of the N complex random variables T^ 
is called trace invariant if its cumulants are linear combinations of trace invariant 
operators, 

C{B) 

K2k [T Hl ,t &i ...T Hk ,f h } = J2 -W^v) II S nn , (3) 

B, k(B)=k p=l 

where the sum runs over all the D-colored graphs B with 2k vertices. 

There exists a unique D-colored graph with 2 vertices (it has D edges connecting 
all the two vertices), called the -D-dipole and denoted B^ 2 >. We are interested in the 
large N behavior of a trace invariant probability measure. In order for such a limit to 
exist, the cumulants must scale with N. We denote the rescaled cumulants 

*&,»*) -K(B,N), (4) 



]\f-2(D-l)k(B)+D-C(B) 

and we call K(B^ 2 ',N) the covariance of the distribution. 

Definition 3. We say that the trace invariant probability distribution is properly 
uniformly bounded at large N if 

lim K(B( 2 \N) = K(B^) < oo , 

7V->oo 

K (B, N) < K(B) , \/B / B {2) and N large enough . (5) 

The scaling in equations ([4]) and ([5]) is the only scaling leading to a large N limit. 
This is a nontrivial statement and the reader should consult [39] for detailed expla- 
nations on this point. The simplest example of a probability distribution for random 
tensors is the normalized Gaussian distribution of covariance a 2 

T D-1 



a 2 2m 



^2 ^— . ( 6 ) 



characterized by the expectations of the connected trace invariants 

Tr 6 (T,f)) CT2 = / (II^^ 1 ) e- ND -^^* T * S ^ Tr B (T,f) . (7) 



The moments of the Gaussian distribution are rather non trivial and have been studied 
in |49j . For any connected graph B with 2k(JS) vertices there exist two non-negative 
integers, 0(B) and R{B) such that 

lim A^ 1+ ^ B )/Tr B (T,f)\ = a 2k ^R{B). (8) 

N^roo \ I a 2 

The normalization in eq.© is the only normalization which ensures that the conver- 
gence order is positive and, more importantly, for all B, there exists an infinite family 
of invariants (graphs B') such that Q(B) = Vl(B'). Again, this is a nontrivial statement 
and the reader should consult H9] for more details. 



Definition 4. A random tensor T distributed with the probability measure fijy con- 
verges in distribution to the distributional limit of a Gaussian tensor model of co- 
variance a 2 if, for any connected trace invariant B, 

lim N^ 1+n ^fi N \Tv B {T,f)] = a 2k ^R(B). (9) 

The universality theorem for tensor models is [39] 

Theorem 1 (Universality). Let N D random variables T^j whose joint distribution is 
trace invariant and properly uniformly bounded of covariance K(B^ 2 \N). Then, in 
the large N limit, the tensor T^ converges in distribution to a Gaussian tensor of 
covariance K(B^) = lim^oc K(B^ 2 \ N). 

2.2 Borel Summability 

The perturbative expansion of tensor models (and of quantum field theories in general) 
is not summable. The root of the problem is that one usually performs an expansion 
in some coupling constant A around A = 0. However the interaction is stable for A > 
but unstable for A < 0. The partition function and the cumulants are analytic in some 
domain in the complex plane outside the negative real axis. Hence A = belongs 
to the boundary of the analyticity domain of the cumulants. A Taylor expansion 
around a point belonging to the boundary of analyticity domain of some function is 
not absolutely convergent. However, in some cases, such Taylor expansions turn out 
to be Borel summable. 

Theorem 2 (Nevanlinna-Sokal, [58]). A function f(X,N) with A G C and N G R + is 
said to be Borel summable in A uniformly in N if 

• f(X,N) is analytic in a disk SKA -1 > R~ l with R G ]R + independent of N. 

• f(X,N) admits a Taylor expansion at the origin 

r-1 

f(\,N) = Y,fN,k\ k + R N ,rW , \Rn,tW\ < Ka r r\\X\ r , (10) 

fc=0 

for some constants K and a independent of N. 

If /(A, N) is Borel summable in A uniformly in N then B(t, N) = J^fcLo vJ N ^ ^ s 
an analytic function for \t\ < a~ l which admits an analytic continuation in the strip 



{z\ |9te| < a 1 } such that \B(t,N)\ < Be l ' R for some constant B independent of N 
and f(X,N) is represented by the absolutely convergent integral 



f(\N) 



dtB(t,N)e~ 



;n) 



That is the Taylor expansion of /(A, N) at the origin is Borel summable, and /(A, N) 
is its Borel sum. The important thing about Borel summability is that it provides a 
uniqueness criterion: if a divergent series is the Taylor expansion of a Borel summable 
function f(X,N) at A = 0, then /(A, N) is the unique Borel summable function whose 
Taylor series is the original series. 

The set {A|KA _1 > R ,R G K+} is a disk (which we call a Borel disk) in the 
complex plane centered at -^ and of radius -^ (hence tangent to the imaginary axis, see 
Figure ED as, denoting A = ■§ + ae* 7 , 



RX 



P 2 

i > fl-i & — > a 



(12) 




Figure 2: A Borel disk. 



2.3 Decorated Trees 

All our results rely on expansions indexed by various kinds of trees which we describe 
at length in this section. We will denote a a -D-tuple of permutations over k elements, 
a = {ai,...a D ). 

Unrooted plane trees with colored, oriented edges and marked vertices 7^f t . 
An unrooted plane tree is a tree with a cyclic ordering (say clockwise) of the edges at 
every vertex. We denote the total number of vertices of the tree by n and label them 
1, 2, ... n. The edges (i,j) of the tree are oriented either from i to j or from j to i 
and have a color c £ {1, 2 . . . D}. Plane trees with marked vertices t = {i\, . . . i^} are 
obtained by selecting a preferred starting point of the cyclic ordering at the vertices 
ii,...ik- The starting point is represented as a mark (or cilium) on the verte^o We 
denote such a tree T^ L with i = {i x , . . . i^} and we denote the abstract tree associated 
to Tni by T n . Note that several plane trees are associated to the same abstract tree, 
and a vertex can have at most one cilium. An example is presented in Figure [3j 



5 Such trees are sometimes called ciliated plane trees. 







Figure 3: A ciliated plane tree with colored edges. 



All the notions we introduce in the sequel refer to a fixed plane tree T^ L - To simplify 
notations we consider this dependence implicit. 

Every plane tree has an associated contour walk (Harris walk or Dyck path) con- 
sisting in the ordered list of vertices encountered when walking clockwise around the 
tree starting from some position. We identify a tree with its contour walk. Whenever 
the walk encounters a cilium, the walk steps twice (first a step to the cilium and then a 
step from the cilium) . We signal the presence of a cilium by a semicolon (rather than a 
comma) separating the two entries in the walk. The contour walk of the tree in Figure 
[3] is therefore 



(1, 2, 3, 4, 3, 5; 5, 6, 5, 7, 5, 3, 8, 3; 3, 2, 9, 10, 9, 11, 9; 9, 2, 12, 13, 12, 14, 12; 
12,2,15,2,16,2) . 

The contour walk of 1~^ L has In— 2+k steps, that is entries, denoted q = 1, 2, . . . 2n— 2+ 
k. The walk in the previous example has 34 steps. We denote the vertex corresponding 
to the step q by i(q). The contour walk is cyclic. 

For every tree T^ L we consider its connected subgraphs having all the vertices of T^ L 
but only the edges of color c of T^ b (see Figure [3|) . Each such subgraph is comprised of 
several plane trees which we denote f c . We call them faces of color c of T^ r The walk 
of T'n L induces a contour walk for each of its faces f c obtained by selecting only the 
(ordered set of) steps q(f c ) at which the contour walk of T^ L encounters a vertex of 
f c . For instance, the face of color 1 with vertices 1, 2, 15 has an induced contour walk 
f 1 = (1,2,2,2,2,15,2,2) consisting in the steps q{f l ) = {1,2,16,23,30,31,32,34}. 
Note that the walk of T^ L can leave and re-intersect several times f c (in the previous 
example it leaves f 1 after the step 2 and re-intersects it at the step 16). We partition 
the faces into two subsets: 

• The faces f c having no cilium. We call them the internal faces of color c of 

T^ L . For the tree in figure El f 1 = (1, 2, 2, 2, 2, 15, 2, 2) or f l = (4) are examples 
of internal faces of color 1. 



• The faces f c having at least a cilium. We call them the external faces of color c 
of Tn L - Say the vertices i^ ■ ■ ■ ik d are the ciliated vertices of f c properly ordered, 
that is f c = (. . . i kl ; iki ■ ■ ■ h 2 ; ik 2 ■ ■ ■ ik d ; h d ■ ■ ■ ) ■ We can further subdivide f c 
into the walks f^H^ik^ f,i k ^i H up to f****^ We call these walks the 
external strands of 7Jp t . For the example of [3l the external face of color 1 
f 1 = (3, 3, 5; 5,6, 5, 7, 5, 3,3; 3) subdivides into the external strands of color 1 
/i;3->5 = (3,3,3,5) and / 1;5 ^ 3 = (5,6,5,7,5,3,3). The strands are comprised 
of the (ordered) entries separated by semicolons in the walks of the faces. Every 
strand f c '< 1 ^ 1 has an induced contour walk, made of the steps q(f c ' l ~ > ' 1 )■ 

Any vertex belongs to exactly D faces, either internal or external, one for each color. 
The cilia are identified as the encounters of semicolons (separating the same label) in 
the contour walk of the tree. We denote qi ■ ■ ■ qk the positions of the cilia in the contour 
walk of the tree 7^° t (that is i{qi) precedes a semicolon and i(qi) = i(qi + 1)). Thus q\ is 
the step at which the first cilium is encountered, q2 the step at which the second cilium 
is encountered an so on. All the faces f c to which the ciliated vertices i(qi) belong are 
external. 

The external faces and strands of color c can be encoded by a permutation £ c over 
1, . . . k. The encoding goes as follows. Consider the ciliated vertex i(q\). It belongs to 
the external face f c of color c. Starting from i(qi), the first ciliated vertex (which is of 
course unique), we encounter in the contour walk of f c is of the form i(qi) for some I (as 
it is one of the ciliated vertices of T^J . We set £ c (l) = I. We repeat the procedure for 
92 and so on and obtain a permutation £ c over 1, . . . h. For the example in figure [3] we 
have q x = 6,q 2 = 14, q 3 = 21, g 4 = 28 (and i(q t ) = 5,i(q 2 ) = 3,z(g 3 ) = 9,z(g 4 ) = 12). 
The permutations £ c write in cycle decomposition as £i = (1,2)(3)(4), £2 = (1)(2,4)(3) 
and £3 = (1)(2)(3)(4). The external faces of color c correspond to the cycles of £ c 



f = (i(qi) ■ ■ ■ *(«f c (o); *(?€c(q) • • • *teo); 

*(ag(o) • • • ifag- 1 ®)'' ifag- 1 ®) ■ ■ ■ *(»)) i Cc(0 = l » ( 13 ) 

and the strands write as 

f i{qi) ^ i{q ^ ) =^(qi)...i(q U l)))> 
f^^^^ffeffll-feio)),.. (14) 

Consider for instance the example of [3] and £1 = (1,2)(3)(4). The external faces and 
strands of color 1 are 

f 1 = (3, 3, 5; 5, 6, 5, 7, 5, 3, 3; 3) => / 1;3 ^ 5 = (3, 3, 3, 5) , / 1;5 ^ 3 = (5, 6, 5, 7, 5, 3, 3) 

f l = (9, 9, 11, 9; 9) =* / 1;9 ^ 9 = (9, 9, 9, 11, 9) 

f 1 = (12, 12, 12; 12) =>- Z 1 ; 12 ^ 12 = (12, 12, 12, 12) . (15) 

We denote the set of all internal faces of T^i by -7 rmt ('7^' i J, the set of all external 
strands of T^\ by 5 ext (7pJ, and the set of all faces of T^ L by F{Tn L )- 

A tree can be built by adding one by one edges connecting univalent vertices. At 
every step of this procedure the total number of faces of a tree increases by D — 1 (every 



^Where we use the cyclicity of the walk to rearrange the last list in the appropriate order. 



vertex connected by an edge of color c will bring a new face for every color d ^ c) , and 
taking into account that the tree with a unique vertex has D faces, we obtain 



D 



\HrZ)\ =D + (n- 1)(D - 1) = |J- int (7^)| + Y, C&) 



(16) 



c=l 



where C(£ c ) denotes the number of cycles of the permutation £ c 



Plane trees 'with external edges. A notion closely related to ciliated plane trees is 
the one of trees with external edges. For every c, we add to the tree dashed oriented 
edges of color c connecting pairs of cilia such that every cilium has exactly an incoming 
and an outgoing dashed edge (we allow an edge to be incoming and outgoing on the 
same cilium). We call the dashed edges external edges. We presented in figured] 
and example of tree decorated by external edges (we represented the external edges of 
color c decorating the faces of color c) . 




<*•/• 



->,■ 





V 




Figure 4: A ciliated plane tree with external edges. 



The external edges of color c can be encoded by a permutation r c over k elements. 
The outgoing external edge of color c on the cilium on the vertex i(qi) is an incoming 
external edge on some other cilium, say the one on the vertex i(qi). We set r c (l) = I, 
and repeat the procedure for all the cilia. 

The example in figure [3] corresponds to the permutations with cycle decomposition 
n = (2)(1,3,4), r 2 = (2)(1,4,3) and r 3 = (1,2)(3,4) (recall that q x = 6,q 2 = 14, q 3 = 
21,(74 = 28 and %i) = 5,i(q 2 ) = 3,i(q 3 ) = 9,i(g 4 ) = 12). 

We denote the tree decorated by external edges 7^^, with r = (ri, . . . m)- The 
external edges recombine the external strands of Tjr, into external faces of T =-. 
They are defined as follows. 

Starting from the external faces f c of T^ L , 

f = (*(«) • • • »(?&(!)); *(?&(0) • • • *(?£(q); 

*(9g(o) • • • ^^-^z)); '(ff^- 1 ©) • • • *(«)) ' ^(0 = z > ( 17 ) 
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with strands 

f#(«Mitcm) = (i(q l )...i(q m j), (18) 

we build the external faces of TP =, indexed by the cycles of r c £ c 

f c ( T Z,r) = (*(«) • • • *'(%c(o); *(9r c &(o) • • • *(flWc(o); 

*(9Tcf c rc€ c (0) • • • *(9fe(r e e e )^i(I))) » (Tc£ c ) d (l) = I ■ (19) 

Consider for instance the example of[U Recall that we have £1 = (1,2)(3)(4), ri = 
(2)(1,3,4) and 

/i;3->5 = (3, 3, 3, 5) , / 1;5 ^ 3 = (5, 6, 5, 7, 5, 3, 3) , 

jl;9^9 = ^ 9) g^ n> g) ^ ^1;12-H2 = (^ 12) 1 2 , 12) . (20) 

We have n£i = (1,2,3,4) and the associated external face of T ? is 

f C (T^) = (5, 6, 5, 7, 5, 3, 3; 3, 3, 3, 5; 9, 9, 9, 11, 9; 12, 12, 12, 12) (21) 

The internal faces oiT nL f are the internal faces of T n ,i and we denote J-(T n L ?) the set 
of all the faces, internal and external, of T n t ?. We have 

c 
D 

= D + (n-l)(D-l)-J2 C(Q + Yl C (Tctc) ■ (22) 



c=l 



Plane trees with external edges and loop edges. The last type of decorated 
trees we will use is trees with loop edges. We add to the plane tree with external edges 
T^ L if, 2s new cilia located on the vertices j\ , jy , ji , J2' up to j s , j s t . The new cilia are 
allowed to be located anywhere on the tree, including on one of the vertices in i. We 
assign to these new cilia a color such that j\ and jy have both color c\ , ji and jV have 
both color C2 an do on. We connect the cilia on j\ and jy by a dashed edge of color c\ 
and so on. We call these new edges loop edges. The loop edges are not oriented. We 
denote the graph thus obtained T^ rc with C = {(jijy), (J2J2'), ■ ■ ■ (js,js')}- We 
present in figure [5l in the upper left corner, a plane tree with two loop edges: one of 
color c\ = 3 and one of color c% = 1 with j\ = 2,jy = 2,J2 = 2,j 2 ' = 3. The loop edges 
have no orientation. 

The presence of the loop edge has several consequences. First, the walk around the 
tree is modified by the presence of the new cilia. For the example in figure [5j in the 
presence of the loop edge, the walk writes 

(1, 20202, 3, 4, 303, 5; 5, 6, 5, 7, 5, 3, 8, 3; 3, 2, 9, 10, 9, 11, 9; 

9, 2, 12, 13, 12, 14, 12; 12, 202, 15, 2, 16, 2) . (23) 

where we boxed the semicolons representing the new cilia (a new cilium transforms 
i, into i;i in the walk). Furthermore the presence of a loop edge modifies the faces. 
Consider the cilia on j p and j p i . They can either 
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V 




Figure 5: A ciliated plane tree with external edges and loop edges. 



belong to two distinct faces of color c p , 
f\ = {■■■3p\3p---) an d 



fCp 

J2 



■ 3p' j Jp' 



(24) 



Then the two faces are merged into the face f Cp = (. . . j p ; j p > . . . j p i ;j p ...). 

• belong to a unique face of color c p . Then f Cp = (. . . j p ; j p . . . j p i \j p >...) splits into 
two faces (j p . . . j p >) and (j p > . . . j p ). 

More precisely, the tree decorated by loop edges T^ L ^ r is the tree with external 
edges 7^p/ =, obtained by adding the 2s new cilia on ji,Jx>,J2,J2', ■ ■ -js,js' (hence J = 
i U {ji , j[ , j 2 , j' 2 , • • • js ,3s}) and connecting by external edges 

• j p with itself (and j p > with itself) by a pair of external edges of opposite orienta- 
tions for every color c/ c p . 

• j p with j p i by two external edges (one for each orientation) of color c p . 

The permutation £' and r' associated to T (J , ™ can easily be identified: for the 
example of figure [5] we have 



q x = 2 , q 2 = 3 , q 3 = 7 , <? 4 
ei = (1,2, 8) (3, 4, 5) (6) (7) 
^ = (1,2,3,5,7,8)(4)(6) 
^ = (1,2,6,8)(3,5)(4)(7) 
r{ = (l)(2,3)(4,6,7)(5)(8) 
r^ = (l)(2)(3)(4,7,6)(5)(8) 
r^ = (l,8)(2)(3)(4,5)(6,7). 



9 , 95 = 17 , <?e = 24 , q 7 = 31 , q 8 = 33 



(25) 



However, writing them requires a bit of care as the presence of the new cilia shifts the 
steps in the contour walk. It is best to present T ,/ » in two stages. The cilia of T! ,,, are 
either cilia of T^\ or they are among the 2s new cilia. The step qi corresponding to the 
Z'th cilium encountered in the walk around T^ L becomes the step q' m n\ corresponding to 
the m{iy th. cilium in the walk around T^\,- The other steps, denoted q' t , ■> and q' t i p i\ in 
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the walk around Tp ' / corresponding to the new cilia j p and j p i . It follows that the steps 
q' r corresponding to cilia in the walk around TP.i are partitioned into three categories: 
r = m{l) for some I, r = t(p) for some p or r = t(p') for some p' . 
For the example of Figure [5] we have 

m(l) = 4 , m(2) = 5 , m(3) = 6 , m(4) = 7 , 

t(l) = 1 , t(l') = 8 , t(2) = 2 , t(2') = 3 . (26) 

Step 1. The presence of new cilia modifies the cycle decomposition of £' c . A cycle 
of £ c of the form (. . . Z£ C (Z) . . . ) becomes (. . . m(l) . . . t(p) . . . t(q') . . . £ c (?n(Z)) . . . ), where 
tip) and t(q') denote the (ordered list of) new cilia inserted on the face corresponding 
to the cycle of £ c (for example the cycle (3, 4, 5) in £[ which comes from the cycle (1, 2) 
in £i). Note that if j p and j q i are inserted on some internal face of 7^p t then they will 
form a cycle . . . t(p) . . . t(q') ... in the permutation £' c having no corespondent in £ c (for 
example the cycle (1, 2, 8) in ££). We reconnect the new cilia with trivial external edges 



T c (r) 



m ( r c(0) if r = m(l) 

t{p) if r = t(p) . (27) 

tip') \ir = tip') 



We thus obtain an intermediate tree with external edges T° — The permutations f c 
acquire extra cycles of length 1 with respect to r c . As a function of the cycle structure 
of £g, one can at most convert internal faces of color c into external faces of colors c, 
(if £' c has more cycles than £ c ) but the total number of faces does not change 



l^fe)! = \^ T Z,f)\ ■ (28) 

Step 2. We convert the permutation f c into the permutations r' c defined as 



T 'c( r ) 



rn ( r c(0) if r = m(l) for some I (t{(4) = 6 ) 

tip') if r = t(p>) for some p and c = c p (r{(2) = 3 ) 

tip) if r = tip') for some p' and c = c p (t{(3) = 2 ) . (29) 

tip) if ?" = tip) for some p and c ^ c p (r^(2) = 2 ) 

^tipl) if r = tip') for some pi and c ^ c p (t^(3) = 3 ) 



We now obtain the tree with external edges 7?, « corresponding to the tree with 
loop edges T^\ - c ■ For every loop edge of color c, two cycles of length one in f c are 
merged into a cycle of length two of t' c . 

On the other hand, for any permutation £, if t' is obtained from f by merging two 
cycles of length 1 into a cycle of length two, 

Cir'O < CifO + 1 , (30) 

hence 

E c «0 < £ c{u' c ) + - =► 1^(7*^1 < |J(^)| + - 

c c 

D 

=* I^CC,^)I < i> + (n - !)(£> - 1) - j; C (^) + E C( ~ T ^ + s ' ^ 31 ) 



c=l 
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where J r (T^ ) L ^ c ) is the set of all the faces of T n ,i,,r,c- 

We will associate to each step q in the contour walk of a tree (or of a tree with 
external edges, or of a tree with external and loop edges) a positive real parameter a q . 

2.4 Interpolated Gaussian Measure 

Consider an abstract tree T n with n vertices labeled 1, 2, ... n. To every vertex 1, 2, ... n 
we associate D matrices (one for each color) of size N x N. We denote the matrices 
associated to the vertex i by a" 1 ,? i • • • c'* • We associate to every edge of the 
tree (i,j) £ T n a real variable u lJ . To every couple of vertices k and I we associate the 
function 

w kk (T n ,u) = l w kl (T n ,u) = inf v*J, (32) 

with Pk^i{T n ) the unique path in the tree T n joining the vertices k and I. We denote 
f^wU (T n , u )i® D ( a ) the normalized Gaussian measure of covariance 

d^ w ij(r n ,u)i^((j) <J i ab C( y (T(l)c ' ] )vai 
= J dn wij(TnjU)l ® D (a) a^ c a%' = w kl (T n ,u) 5 aa ,5 w 5 cc ' . (33) 

The existence and uniqueness of this measure follows form the positivity of the real 
symmetric matrix w (T n ,u) (see eg [59]). 

2.5 The Graph of D Permutations Over k Elements 

Any D-tuple a can be represented as a D colored graph B# with labelled vertices. 
We draw k black vertices labeled l,...k and k white vertices labeled 1, . . . k and we 
connect the vertex I to the vertex a c {l) by an edge of color c oriented from I to cr c (l), see 
Figure[6l Conversely, to every colored graph B with 2k vertices labelled 1, ... k, I, ... k 




*c(l) <r c (l) 



Figure 6: Graph associated to a D-tuple of permutations. 



we associate a unique D-tuple of permutations <?(<6) encoding the connectivity of its 
edges. 

2.6 The Weingarten Function 

The Weingarten function introduced in [54^ [55] arises naturally when one considers 
integrals over the unitary group, namely 

/ [dU}f[U n]P] Ul nl = 

JU{N) fj[ P > i 
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= E <W ctC1) • • • S nuK (k) S plP ' Tm ■ ■ ■ S pkP ' T(k Wg(iV, ar 1 ), (34) 

where the sum runs over all the permutations a and t oik elements and the Weingarten 
function is 

where the sum runs over the partitions tt of N, x n is the character of the symmetric 
group corresponding to n and s 7r) Ar(x) is the Schur function of the unitary group (hence 
s 7r,Ar(l) i s the dimension of the irreducible representation of U(N) associated with tt). 
We will in particular use the following properties of the Weingarten function [Ml E5] 

Wg(iV,(l)) = l, 

lim N^~^W g (N,a) = ff (-l)l*.MI-i * ( Wf)} " 2 ) , (36 ) 
Af^oo ax |C s (cr)| V |C s (cr)| - 1 / 

where C(<r) denotes the number of cycles of the permutation a and |C s (cr)| denotes the 
length of the s'th cycle. It follows that for N large enough we have 

IWg^^K^cH 22 "' ( 37 ) 

and we will always assume in this paper that N is large enough such that this bound 
is respected. 



3 The Quartically Perturbed Gaussian Measure 

Our starting point is the quartically perturbed Gaussian tensor measure 

^(4) = ^^(n^ 1 ^^) e- ND - lsW V,f) t (38 ) 



D 

\ "* T T T T X X TT X X 



5 (4) (T) f ) = £ TnS n n% + A E E T A T A W*m<»< II <W^ 



with Z(A, AT) some normalization constant. The quartic perturbation corresponds to 
a sum of invariants whose graphs are represented in figure \7\ 



D ••• 2 2 (■•■ D 



Figure 7: The graphs of the quartic perturbation terms. 



One can in principle consider more general quartic perturbations (in the associated 
graphs the vertices share q and D — q lines respectively with q > 1) however such 
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perturbations are suppressed in powers of 1/N. The generating function of the moments 
of /jt 4 ' is 

Z(J,J;X,N) = [ {J\N D - 1< ^^-) e -^- 1 ^>(T,f) + E,f,j, + E rI T„J„ ; (39) 

n 

and the generating function of the cumulants is W( J, J; A, N) = In Z(J, J; A, N), 

Q(2k) 

«(T*.T fc ,...T« b ,T Sl ) - w - BJ -^ J - JiJ -W(J,J-,X,N ) \ j ^ o . (40) 

Our first result concerns the constructive expansion of W(J, J; A, N). 

Theorem 3 (Constructive Expansion 1). The generating function of the cumulants of 
/i( 4 ) in eg. (J38j) writes as a sum over plane trees T^ t 

1 n 1 1 n 

W(J, J; A, N) = j: ^(-A)- 1 E fc| jyO^D^i) E E^O'W 

n>l fc=0 H^2--'. l fc= 1 T„ t 

where the contribution of a tree is 

2n-2+k 



l Z{T%J=J ( II dut3 ) J d ^{T n ,u)i^(a) J [ J] <W)e~ E "" 



(Jj)ST n " " q=l 

/c e jrint (Tn Oj g 6 t/(/ c ) {pf,nf} 1=1 

r ( c r i (9)c_ (T i(9)ct)- 



X 



n n 



where we used the notations of section \2.3\ and — > means that the products are ordered. 

This expansion is the generalization to tensor models of the constructive Loop 
Vertex Expansion (LVE) introduced in |51j for matrix models. Note that our expansion 
looks somewhat different from the LVE of [51] , notably the step parameters a p have no 
equivalent in the initial formulation. The usual LVE is recovered from [3] by restricting 
to trees having no ciliated vertices. In this case the integrals over a p can be computed 
explicitly and one recovers the formulation in terms of resolvents of [51] . 

As a consequence of the (LVE) we can derive an expansion for the cumulants of our 
measure. 

Theorem 4 (Constructive Expansion 2). The cumulants of the measure fi^' in eq. ([38]) 
are trace invariants 

C{B) 

«(T ft ,fa 1 ,...T A ,f«J= £ £(*?,/$) !]<*£> (42) 

B, k(B)=k p=l 



1(3 



where the sum runs over D-colored graphs with 2k vertices labelled 1, ... k, 1, ... k and 
admit an expansion as a sum over plane trees with external edges T n L ? 

1 i n 

n>k i 1 ,i 2 .-,i k =l'j<> l ? 

where the contribution of a tree with external edges is 

D /"l /■ 

'Z £ (Tn ) ,J = kl(j[Wg{N,T c a- 1 (B))) / ( J] du^) d^ (TntU)1 ® D (a) 

1 J / - _-\ _ m "^ 



c=i ° (i,i)eT n 



2n-2+fc 



x/( n ^>- E ^ 2+fea * n t 4 n e - Q '^> i(,w(9)ct) 

The main advantage of the constructive expansion of the cumulants is that, unlike 
the perturbative expansion in A, it leads to a series which is absolutely convergent 
uniformly in N. 

Theorem 5 (Absolute Convergence). The series in eq. (|43p is absolutely convergent 
for A 6 1, A G [0, 2~ 3 D~ 1 ). Moreover, the cumulants are bounded by 

\&(B,$)\ < N^^^^^lX^KiB) , (44) 

for some constant K(B) independent of N (and independent of X for |A| small enough). 

The scaling bound of equation (I44D coincides with the proper uniform bound of 
definition^ However we can not yet conclude that fi^ A ' is properly uniformly bounded. 
Indeed, in order to conclude this, one still needs to prove that the second cumulant 
converges when iV — > oo. We will show this later in this paper. 

The domain of convergence < A < 2 _3 L> _1 is optimal. We will see below that the 
leading order in 1/N, the "melonic" sector yields a series whose radius of convergence is 
exactly 2~ 3 D -1 . It is not surprising that the full non perturbative expression diverges 
when one reaches the radius of convergence of the leading order. 

The series (|43|) computes the cumulants for real, positive (and small) coupling con- 
stant. The cumulants can be analytically continued to some domain in the complex 
plane. 

Corollary 1. The cumulants K(B,fi N ) can be analytically continued for A = |A|e lI,p 

with ip 6 (— 7T, it) and |A| < (cos ^) 2~ 3 D~ 1 . In this domain they are represented by 
the absolutely convergent series 

m^ ) ) = T / ^r- 1 w ^ m z T) t ££^(0> ( 45 ) 

D 



^(T^) = k\{j\Wg{N,r c a-\B))) !\ J] duv) [ d» wij{Tn>u)1 ®n(a) 

c=l '■ *~™ 

/2n-2+k 
( e~ l ^da.q\e~ 



c=l J ° 

2n-2+k 

HI 

9=1 
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x H Tr[ f[ e- a ^T& 



j^^Wc-^Met) 



We have thus a well defined expression for the cumulants in a heart-shaped do- 
main (see Figure [8] below where we represented by a dashed circle the circle of radius 
2~ 3 L>~ 1 ). As expected, A = is a point belonging to the boundary of this analyticity 
domain. 




Figure 8: Domain of convergence. 



The interplay between the non perturbative LVE expansions presented so far and 
the 1/N expansion is captured by the following theorem. Consider the rescaled cumu- 
lants, which according to eq. (j4"3j) write as 



,(4), 



K (B,N) = N- D+2k ( D -V+ c W&(B,fx k ' 



(46) 



n>k 



'i' i 2-.')t= 1 75J l ? 



where the rescaled contribution of each tree with external edges T - is 

T E (T^J = N -D+2k(D-l) + C(B)-(k+n-l)(D-l) fc , (j[ W g(iV, T c O-\B) 

c=l 

x / ( II dn ' J ) d ^w»(T n ,u)i^{a) 



(»J')6T„ 

2n-2+fc 

n 



^2n-2 + fc -i-* 



e^ e -E^^«. 



x n t 4 n e -^v^(- i(9)c - iwcf ) 



(47) 



Theorem 6 (The Mixed Expansion). 27ie contribution of a tree with external edges 
Tn L admits an expansion in terms of trees with external edges and loop edges T^ > L - c 

9=0 



t (m) ^J= E Y. TELM Vn M ) ( 48 ) 

C,\C\=q ci...c g 

where C runs over all possible ways to decorate T n f * with q loop edges (ji,ji'), (J2,J2') 
U V to (jq,j g i) and c\,...c q run over the possible colorings of the loop edges and the 
contribution of a tree with external edges and loop edges is 



c=l 

-1 s 

C I ( du lJ ) J 



n d^)n^v(r n/ 



(tj)eT„ p=i 

xiV -D+2fe(D-l)+C(B)-(fe+n-l)(D-l)+|^(T„» i - i£ )| ^ ^ 



while the rest term is 



R {s \rZ,r)= / dt{i-t) 

JO 

1 / A 



s-l 



N -D+2k(D-l)+C(B)-{k+n-l)(D-l) 
D 



)IV ^'Kn^WiB))) E E 

V ; ' c=l £,|£|=s ci...c s 

' ( n d« y )n«^^(rn,«) / dA*wii(T n ,«)i^(ff) 

(ij)eT„ p=l 

oo 2n-2+fc+2s v 

nip \ \-^2n — 24-k-\~2s — i-?f 



x 

3=1 



n Tr f n e - a ^^^T (<TiC9)c_(ri(9)ct) 



(50) 



Furthermore the terms in the mixed expansion admit the bounds 

IrfcJfT il < |A| " (/: e° fc W 2n + 2g + fe ~ 3)! f5i) 



Ni( D ~ 2 ) y ' y g!(2n + A;-3)! 

1 V n ^ T>] ~ ( cos |^n+2,+fc-2 #,(0-2) <> • >( s _ 1 )l( 2n + jfe _ 3)| * 

We call this expansion the mixed expansion because it is at the same time an 
expansion in A and an expansion in 1/N. More precisely, being an expansion in N u-i , 
one can use it to establish the Borel summability of the cumulants or, alternatively, 
one can use it to establish the -^ expansion of the cumulants at all orders. 

Theorem 7 (Borel Summability). The rescaled cumulants 

K(B, N) = N -D+2k { D-l) + C{B)^ B ^ ^ ^ (52) 

are Borel summable in A uniformly in N. 

A crucial point is that, as we are interested in the large N regime, both the conver- 
gence of the constructive expansion in its analyticity domain and the Borel summability 
around A = are uniform in N. 
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Theorem 8 (The 1/N expansion of the cumulants). Using the mixed expansion, the 
rescaled cumulants of /j,^' write as 

K ^ N ) = E^r 1 E EEE E Y, tELM (tZ,?,c) 

n>k i i> i 2---. i fe= 1 T„° t ? q=0 C,\C\=q ci...c, 

+R { x\B,\), (53) 

and for |A| < 5 -2 2 -1 .D _1 (cos tj) t/ie res£ £erm admits the bound 

\R%'(B,X)\ < Ka S S \ jj—rr ' ' - J _„. : . (54) 

/or some constants K and a. 

One can now use the 1/iV expansion as follows: the terms up to order N~ S ( D ~ 2 ' are 
indexed by trees with at most s loop edges. Such corrections can be evaluated order 
by order. In particular all corrections at fixed order in 1/N will reach criticality when 
the sum over n becomes critical, i.e. all terms in the 1/N expansion will diverge for 
the same critical constant. 

The factorial bound in eq. (|54p suggests that the cumulants are Borel summable in 
1/N. This is most likely the case, however the attentive reader will notice that we 
did not yet establish analyticity of the cumulants in 1/N. This is difficult, because, 
besides the explicit occurrences of N, we also must take into account that N is the size 
of the matrices ac ■ In order to establish analyticity in 1/N one needs to find a better 
representation of the cumulants in which TV appears exclusively as a parameter. We 
postpone this to future work. 

Before we conclude this paper we present as an example the leading order behaviour 
of the two point cumulant. A number of simplifications arise in this case : k = 1, 
Wg(iV, (1)) = -^ , C{B,^S) = 1. Using the 1/N expansion up to order s = 1 we get 

1 n 

^ (2) ^) = E^(-^- 1 E E tEL ' {0) K {hU ^^ r n^ 2) ^^ 

n>k 11=1-7-0 

'n,{i x } 

t el,(o )(t o ) = J_ N -D+2(D-D + i-(i+n-i)(D-i)+\m^ ^ 9 )\ 

{ n,{nUlW N D 

iV ( cos f ) 

The leading order behavior can be resummed, as 

LF(T° - )| = D + (n-l)(D-l)^T EL ^(T° -J = l, (56) 

and the sum over trees can be computed explicitly (we do this for arbitrary trees in 
equation (j73j) below), 

K {B V,N) - YJMW-.M + J,,W . (57) 



n>k 



(n -l)lnl N»-z ( cos |)' 
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Corollary 2 (The large N limit). For - — ' ' 2 small enough 

(cosf) 

^K(BW,N)= - 1 + V1 + 8DX . (58) 

In particular this coupled with the uniform bound in theorem [5] proves that the 
measure /j,( 4 ' is properly uniformly bounded, hence according to theorem [H becomes 
Gaussian in the large N limit. 
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4 Proofs 

Before going to the core of the proofs of the various theorems in the text we establish 
a number of results we will use recurrently in the sequel. 

4.1 Technical prerequisites 

The main inequality on permutations. Recall that C(a) denotes the number of 
cycles of the permutation a and C{Bg) denotes the number of connected components 
of the graphs associated to the D-tuple of permutations a. 

Lemma 1. Let ^ and a and f be three D-tuples of permutations over k elements such 
that £ are permutations encoding the external faces of a tree. We have the bound 

Y C(t c o-^) - Y C{Q + Y, C(t c £ c ) <(D + l)k - C{B S ) . (59) 

c c c 

This bound is a trivial consequence of the following two propositions. 

Proposition 1. Let £ and a be any two fixed permutations over k elements, and let r 
be any permutation over k elements. Then 

CiTo-" 1 ) + C{tO < CiC 1 ^ 1 ) + k . (60) 

Proof: The bound is saturated if r£ is the identity permutation. 

Consider then the sum C(tcf~ 1 ) + C{t^) and suppose that r£ is not the identity 
permutation. The cycle structure of the permutation r^ can be easily read of by 
drawing the graph of the two permutations £ -1 ,r (where we reverse the orientation 
of the edges representing the permutation £ -1 ), see Figure [9j a moment of reflection 
reveals that the number of cycles of r£ is the number of connected components of 
the graph as the t£ jumps from a white vertex to the next white vertex following the 
arrows. 
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r(m) 



Figure 9: Graph associated to £ 1 ,r. 



If r£ is not the identity, then there exists p such that £(p) 7^ r 1 (p) (and r(£(p)) 7^ 
p). We compare C^ff" 1 ) + C(t£,) with CVcr" 1 ) + C^), with t' defined as 



r'(g) = 



r(q) 



Mq^T 1 (p),£(p) 



(P)) = r(e(p)) 



(61) 



U'(£(p))=p 

The graph of the permutations £ _1 ,t' is represented in Figure [H)l By substituting r 

r(£(p)) 




Figure 10: Graph associated to £ , r'. 



with r' we created a new connected component, C(t'£) = C(r£) + 1. Representing now 
the graph of the permutations a,r and <7, r', we see that replacing r by r' amounts 
to permuting the endpoints of the edges representing the r permutation originating in 
T~ l {p) and £(p), hence the number of cycles of rtr" 1 can not decrease by more than 1. 
Thus 



C{t<j~ 1 ) + C(tO < Cir'a" 1 ) + C(t'0 . 



(62) 



Iterating, we find that the maximum is achieved for r such that r(£(p)) = p for all p. 

D 

Proposition 2. Let £ and a be two D-tuples of permutations (with D > 2) such that 
£ are permutations encoding the external faces of a tree. Then 



E c ^) - E c ^) + °( B ^ ^ k - 



(63) 



Proof: For every cilium I we draw a black and white vertex and represent the graph 
Bg associated to the permutation a, see Figure [TTJ 

Consider the tree with associated permutations £. As the univalent vertices with 
no cilia of the tree have no bearing over the permutations £ we can eliminate them. 
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Figure 11: Ciliated vertex and the associated vertices of Eg. 



Consider a ciliated univalent vertex Z in the tree, and say that the line touching it has 
color c. It follows that we are in one of the two cases 




£ c (0 = I , 



6* (0 

Ccx(0 



l Vci 7^ c, 
Z Vci 7^ c. 



(64) 



The permutations a fall in one of the following three categories 



case 1 <r c (l) 7^ I , Vc , 

case 2 a Ci (1) = I , i = 1 . . . q < D , a Cj (l) ^ I 

case 3 cr c (Z) = Z , Vc . 



{*}U{c j } = {l,...D}, (65) 



We build the graph obtained by eliminating the ciliated vertex Z in the tree (by cutting 
the tree line touching it) and deleting the vertices I and Z in Bg and reconnecting the a 
lines respecting the colors (see Figure ITT1 on the right). The new graph is characterized 
by permutations £ and a' over {1, . . . k} \ {1} 



easel e c (p)=Up),Vp^& 1 (l) 
case II £ c (p) = Cc{p) , Vp / Z 

e ci (p)=u(p),vp^i 

case 1 a' c (p) = a c (p) , Vp / cr~ x (Z) , 
case 2 a' c . (p) = a Cj (p) , Vp / cr" 1 (/) , 

<(p) = ^(p) ,Vp/Z 
case 3 c' c (p) = c c (p) , Vp / / 



fi(c 1 (0)=ee(o 






a cA% 



x (0)=^c(0 

1 (0)=^(0 



(66) 



- Vc 



Consider now the graph of two permutations £~ 1 ,o" c . The permutations (££ 
correspond to the graph whose vertices I and I have been deleted and whose edges have 
been reconnected coherently. The number of cycles changes as 

• if a c {l) + I, and &(Q / I then C(a c £ c ) < C(o^ c ) + 1, 
. if o- c (0 + I, and £ C (Z) = I then C(a c £ c ) = C(<7&), 

. if a c (l) = I, and &(Q / / then C(a c £ c ) = CK£)> 

• if tr c (0 = Z and £ C (Z) = Z then C(a c £ c ) = C{a' c Q + 1. 
Finally, the number of connected components of Bg changes as 

• case 1, Bg < Bg/ + 1, 
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• case 2, Bg = Bg/, 

• case 3, Bg = Bgi + 1. 

As a function of the case we are into, we then have 



n 



£C(<r c & 

c 

771: £C(ff«£ 

c 

72: 5^c(o- c e 

c 

7/2: ^C(a c e 

C £c ( £ 

c 

£c(a c e 

E c (e 

c 



13 



7/3 



c 

=X)c(£)+i>-i, c(^)<c(^)+i, 
=£ck&, 

c 

= £C(£) + D, C(B*)<C(B*) + 1, 

c 
c 

= £C(£) + I>-1, C(B 9 ) = C(B 9 ,), 

c 
c 

= £C(£) + I>, C(B 9 ) = C(B9>) t 

c 

= £ C(£) + £> - 1 , C(^) = C(By) + 1 , 

c 

= £c(a^) + 7J, 

c 



(67) 



In all cases 
£ Cfea" 1 ) - £ C(&) + C(B„) < E C^) -1 ) - £ C(0 + C(M + 1 , (68) 

c c c c 

and the bound can be attained only for 73, 773, 72 with q = D — 1, or 71 if D = 2. We 
iterate the procedure, taking into account that after eliminating the univalent vertex / 
of the tree one might need to eliminate some new univalent vertices with no cilia (as 
the latter have no bearing on the permutations £'). Iterating up to permutations of 1 
element £ c = a c = (1) Vc, we get 



£ Cfa- 1 ) - J2 C(&) + C(B a ) < (k - 1) + 1 . 



(69) 
□ 



Evaluating derivatives. We will repeatedly use in the sequel the following result. 
Lemma 2. For any N x TV matrices F and G and any function H(ai,at2) we have 



/•oo 

/ daida 2 H(a 1 ,a 2 ) Y] F niPl G n2P2 

Jo n ,Z~Z„„ 



Pinip 2 n 2 
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Y — 



-ai^/^t^Wc-aWct) 



p 1 n 1 



T 0')=t r 



^-(o-OOc-crOOct) 






p-in-2 



D - / d/3id 7 id/? 2 d7 2 #(/3i + 71 , /3 2 + 72 ) 







xTr 



Fe 



-^V^^^ -^^) -72^^( CT(j)c - CT(j)ct ) 



Ge 



-ft^^(cr«c- CT W)ct) ^yZT^Mc^Mct) 



(70) 



Proof: Expanding the exponentials and evaluating the derivatives the integral writes 

(71) 



poo 

I daida.2 H{ai,a,2) \] F nipi G. 
Jo „,±^,„ 

E 



nil 



E 



pinip 2 "'2 
(- ai ) Ul+Vl+1 



(ui + vi + iy. v a^- 1 



ui : ui>0 

[(ct« c - a« ct ) ui 



™2P2 
«1+«1+1 



Pi" 



r a « c _ a «4w 



611 1 



E 

«2,*>2>0 



(_ a2 )«2+,2 + l / 



(u 2 + v 2 + i)r '\ n d 



-"2+^2+1 



\(a^ - a 



(i)ctl«2 



P2&L 



(o 



U)c _ a U)c\\V2 



an-2 



Using t — , * ppri = f n dx\± — ^V — -V this rewrites as 



N D 



/•oo rl 

— I daida.2 H {0.1,0.2) \ a\dx\02dx2 YJ Fn lPl G n2P2 



E 



pinip 2 «2 

-ai(l - xi)]" 1 [-aixi] Ul [-q 2 (1 - x 2 )] U2 [-a 2 X2] V2 



U\,Vi,U2,V2>0 

— Ul +«1 +«2 +«2 

A 



>'i 



"2 



N D-1 



r (i)c 



a y 



{i)o\\Ul( (j)c _ rr {j) C '{\V2 



(v 



a' 



>'2 



Pin.2 



(T 



U)c _ a (j)cUu2 r a (i)c _ a {i)cjyi 



(72) 



P2H1 



and changing variables to /3i = ai(l — xi),7i = a\x\ and similarly for 2 and summing 
over u\,vi,U2, i> 2 the lemma follows □ 

Combinatorial countings. We count the number of plane trees with n vertices, 
k ciliated vertices i\,...ik and colored oriented edges. Every combinatorial tree T n 
with degrees of the vertices d\, . . . d n , has (2D) n ~ 1 di 1 \ . . . di k \ fli^i (di — 1)! associated 
plane trees with colored oriented edges and marked vertices 1~^ J L , corresponding to the 
two possible orientations of the edges, the D possible colorings of every edge and the 
permutations of all but one of the halfedges touching each vertex (plus a choice d{ r of 
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where to place the cilium on the marked vertices) . The number of combinatorial trees 

i-^Cfl-iy. and we § et 



ith assign 


ed degrees 


d\ , . . . d n 


is 






fixed 

n 






= (2D) n - 


1 E 








d 1 ,...d n = 

J2di=2n- 


i 

-2 



((il -^:.K-i)i d " ! - d "'n'<- i > ! 

= (2Dr 1 (n-2)! £ ^•••^ = (2^r- l ( ( 2 ;+^)! ' (73) 

d 1 ,...d n — l 
Y / d i =2n-2 

as the sums over di yield the coefficient of the term of degree x 2n ~ 2 in the expansion of 

^■l-x' \ (1 - x) n ~ k ~ (1 - x) n+fc ^V P / 

For every tree T^ L one has /c! D trees with external edges T^ L -, 

yi)\ , =k\ D . (75) 

^ / l7# 4 fixed V ^ 

T 

The number of plane trees decorated by loop edges is counted as follows. The walk 
around 7^f t has 2n — 2 + fe steps. The insertion of s loop edges consists in the choice of 
2s positions to insert new cilia. One can insert a new cilium at every step of the walk. 
However, with every insertion of a cilium the walk acquires a new step, hence one has 

[2(n - 1) + k][2(n - 1) + k + 1] . . . [2(n - 1) + k + 2s - 1] = ^^g^' , (76) 

ways to connect the 2s loop edges on T^l- Furthermore for each choice, one has D s 
possible colorings of the loop edges, hence 

E E = D>( 2n f + 28 + k -W . (77) 

CM~.CZ. 'W*- i^n + k- 3)1 

Finally, 

n 

n! ^— ' in — k)\ 

i 1 ,i 2 ---,i k =l 
id+i-d' 

Computing the logarithm of a Gaussian integral. The logarithm of a Gaussian 
integral can be computed using the universal Brydges-Kennedy-Abdesselam-Rivasseau 
forest formula [59] and a replica trick. 

Lemma 3. Let X be a complex vector of components X\, . . .Xn and let I be the 

Gaussian integral of covariance C defined as 

I = j dfic(X) e v ^^ , J duc(X) X a X- b = C al . (79) 
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Then 



ln/ = E^EA II duij ) I '^cr.,»)c(* (0 ) 



n>0 T„ 



(i,j)eT n 



II E ~7Z^) Cal ^U) + E T^m* 7 ^ 



(tj)6T„ \ aft C ' A " ^fe 



ax® ab dx£\ 



~[V(X®,X®) 



(80) 



i=l 



where T n runs over all the combinatorial trees with n vertices labelled 1, 2 up to n, 
(i,j) denotes the tree edge connecting the vertices i and j, the parameters w l: >(T n ,u) 
are defined as 



w ll (T n ,u) = l, w l] (T n ,u)= inf u k 

(*,06P^i(T n ) 



u 



(81) 



where Pj_>j(T n ) denotes the unique path in the tree T n joining i and j and the interpo- 
lated Gaussian measure dfi w ijr Tn U \ C (X^') is 



dW(T n ,u)c(X®) X^X^ = w^(T n , U )C ah 



Proof: We Taylor expand in V(X, X) to get 

/= f dnc{X) e v &>^ = i d^ c {X) Y^~,V{X,X 



(82) 



(83) 



n>0 



The term of degree n can be rewritten as a Gaussian integral over n replicas X^ 1 ' , X^ 2 ' 



up to X' n ) with degenerate covariance between the replicas C r = C ab ; 



ab 



1 /■ n 

1 = E - / dn c ^){X®) \V(X^,X^] 



(84) 



Each term in this expansion is a function of parameters x 13 = x Jt , evaluated for x lJ = 1, 
corresponding to a Gaussian measure with covariance 



^ab - ^ab > 



C 



ab 



x C ab , if=j. 



(85) 



Consider n vertices labeled l,2...n and a function / depending on ^^ — - edge 
variables x lJ with i ^ j. The universal Brydges-Kennedy-Abdesselam-Rivasseau forest 
formula 1591 states that 



/(i,...i)=E/ n «*«*) tt — 



Q\£(Fn)\f 

dx l J 



x kl =w kl (F n ,u) 



,.kl 



W 



(F n ,u) 



inf 



i'j 



(86) 



where F n runs over all the forests built over the n sites, |£(F n )| denotes the number 
of edges in the forest, Vk^l{F n ) is the unique path in F n joining the vertices k and I, 
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and the infimum is set to zero if k and I do not belong to the same tree in the forest. 
In order to apply the BKAR formula to the term of degree n we evaluate 

d\£( F n)\ 



n 



(i,j)e^n 



dx l i 



dwcjx®) n^ w '* w ) 



0|£(*VO 



n 



(ij)e^n 



Q x l] 



E Q 



a b 



tC„ 



tC„ 






Jy(i< ! U (i) ) 



*=i 



E Q 



,6,i (i)^ab -li)~ r '/-,a.b,i^j ■*■ (i)'- / ab ^777 



x«=x«=o 

d 



n (e 



a 



-a 







ab 



ax« ab dx^ ^axP 



+£ 



a a 

^ab 



ab 



ax 



w 



/ 



d^x^c-AX 



(*)< 



]V(x«,x«) 

i=l 

n (e 8 



(i^eJ^n ab 



dx% 



W 



xM,xW=o 



ax 



n+E 



(J) 



5 fl 



(i7< 



axi 



(i) 



ax 



(<) 



Jf(xW,x«) 



(87) 



where we have taken into account that x lJ = x 31 . Thus 

HF„,u)cJX m ) 



^e^e/'( n 
n (e 



n>0 
X 



(M')eF„ 

a a 



(«j')eF n 



ax« 



ax 



CO 



^«(F„,u)C oi 

+ E 



dX a ~ 



d d 

~^ab 



(i) 



ax 



w 



Jy(x»,x 



(»> 



i=l 



The lemma follows by noticing that the Gaussian integral factors over the trees in the 
forest and recalling that the logarithm of a function which is a sum over forests of 
contributions factored over the trees is the sum over trees of the tree contribution. 

□ 
The most important feature of the BKAR formula is that the matrix w 1J (T n ,u) is 
positive |59j . The Gaussian measure is thus well defined and the expectation of any 
function of X, X is bounded by its supremum. 

4.2 Proofs of the Theorems 

In the reminder of this section we present the proofs of the theorems enunciated in the 
text. 

4.2.1 Proof of the first constructive expansion theorem [3] 

The Loop Vertex Expansion of W(J, J; A, X) in eq. (|4ip is obtained by combining the 
Hubbard Stratonovich intermediate field representation and the BKAR formula. We 
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will sometimes drop the bar over the indices of T. 



Step 1: Hubbard Stratonovich intermediate field representation. For any complex num- 
bers Z\ and Z2, e _ZlZ2 can be represented as a Gaussian integral 



dzdz 

2tTT 



e -zz-zZ!+zZ2 z=x+iy f dxdy_ e - x Z-y2- x t y z 1 -Z 2 )-iy{Z 1 +Z 2 ) 



e 4 



z=x—iy I 7T 
-Z1Z2 



(89) 



It follows that a term in the quartic perturbation of the measure in eq. (|38p can be 
represented using N x N integration variables a c ah as 

Tj <6^ ) e -E a6 ^ 6 

aft 

-D-l / \ _ O-l 



-VX2V-T- E(^%n c vc v«c) CT ^=+^ Ar ^~ slwMn^v^) 8 '* 



(90) 



The new integration variables a c ah form D matrices of size NxN, known as intermediate 
fields. Denoting 1 the identity matrix of size N x N we write more compactly 



e ~N D X AEf = l T,nn T R T rh T rh^ii S n c^cS m c n c U c 'jt c S n c' s c'* m c' m c> 



(91) 



f (yi <KtfKb \ g -trK^t ) _^-y_ i A_ Elr ^ Ecl8 c-l 8(ff c_ (T ct )8l8D -cj fs 
caft 



thus Z( J, J; A, N) becomes 

Z(J,J;X,N)= [ (jJjV 



Z) _ 1 (iT^rfT^ 



2vu 



n 



da ab d(7 ab \ p -Ec tr ( CTCCTCt ) 

2va 



(92) 



Xf 



a^- 1 ET« 19B+J-* e c i® c - 1 ®K- f r c t)cg>i® D - c f-+ET-J-+ET s J s 



As cj — <jT is anti hermitian the operator R(o~) 

-1 



1®D _|_ 



7V^ 



-Eci®^ 1 ® (^ C 



c t) (8) 1« 



iD-c 



, which we call the resolvent, is well defined. The intermediate field 
representation renders the integration over TT Gaussian, thus 

Z(J,J;\,N) = /(n d l^°O e " EctoCCTCt+trto(fl( ^ +] ^ <J|R(<T)|J> ' (93) 

c;ab 

where (J\R(a)\J) = Eg,™ ^ fi ( (T )nm4^ 

Step 2: Extracting the logarithm. Using lemma [3] the generating function of the con- 
nected moments W( J, J; A, N) is 

W(J,J;\,N) = J £^J2 I ( II du%3 ) [ d ^ J{ T n ,u)i^(a) (94) 
n>l n ' T n Jo (i,j)eT n 



D 



11 zL ( Z> 8 (i) c 8 jjU + ^ 8rr U)c B Wot 

(ij)eT n c=l V aft ^aft ^fta aft ^aft ^fta 
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n 1 



where the interpolated Gaussian measure writes formally as 

rf Mt U ^(T„,«)l®- D ( cr ) = 



(95) 



Z_*c— 1 I *— 'i t—'a 



a a , a a 

JOc ~(Oct + ~Wc — Wet 



+Z_,j<j x 2^,ab 



a j0c a Ji) C t + SCT U)c s J) C t 
cio oa ab ba 



Expanding the product over i, we get 

W(J,J;\,N) = J2^J2J ( II Mi) J d» wij{TntU)19 v(a^) (96) 



n>l '"" T„ - 70 {i,j)£T n 
D 



11 Z_> Z> Q _« c Q _(i)ct + ^ ' 



(ij)eT n c=l ^ ab 0(J ab °°ba ab a ° ab a(J ba 

n -. -. n 

fc=0 «l ,«2 ■■■^k~ 1 

n 

x J] trln[i?(a (i) )] . 
i^il,...ik 

We need to evaluate the action of the derivative operators on the product resolvents. 
This is done in Lemma [H below. For each combinatorial tree T n one obtains a sum over 
all the plane trees 7* with colored oriented edges compatible with it, and indexing 
the sum by these plane trees we get 



n n 

W(J, J; A, TV) = £ -(-A)- 1 E fcT j^n-ixs-i) E E 

n>l fc=0 i 1 ,i 2 .--,tfc=i 7^\ 



(97) 



2n-2+fc 



/ ( n dt * w ) / d A*to«(Tn,«)i® D ( ff ) / ( n ja 



E2n — 2+fc 
9=1 "8 



n Tr [ n .-^^En^^ 



pn Z=l 



n 



n 



l^-^=T^ iiq)c -^ 9M ) 



^;<(9i)^*(9e c Ci)) 6 5»rt(7^ t ) ?6 g(/ Cii( " ) "' i(? «<=(0 ) ) 
which proves the theorem. 



«■»&(!) 



D 
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Lemma 4. The contribution of a tree T n with marked vertices i±, ■ ■ -ik is 
11 Z> ( 1^ TTWT Q _(.?) c t + ^ ' 



(»,j)6T„ c=l v ab ° a ab ° a ba ab ° a ab °°ba 



(98) 



(J\R(a^)\J) . . . (J\R(a^)\J) trhx[R(<r®)] 



i-xr- 



,. 2n-2+k 

^nyE ( n da « 

T° 9=1 



1=1 



E2n — 2+fc 
q=l a 9 



n Tr [ n .-^■ ,< ' , - , ' M, ]En^..^. i ,.. ? 



pn 1=1 



n 



n 



-<V^T ( ' j!<9)C - CTl(9>Ct) 



/ c;i(9i) ""* (9 5c(0' 6 5^t(^) ge<?(/ c;i(<? ' Mi( ^W ) ) 



P? n £c(0 



where the sum runs over all the plane trees with colored oriented edges and ciliated 
vertices ii,...if~ which reduce to the combinatorial tree T n . 

Proof: We orient the edge (i, j) of the tree from i to j for the term d a (i)d a (j)]. Taking 
into account the sum over c we obtain the sum over trees with colored, oriented edges. 
In order to compute the contribution of each such tree we need to evaluate the action 
of the derivative operators on the product of traces. 

We set Af(ffW) = [U(<7«)]-i = 1® D + J^hr E?=i l® c_1 «>(<x« c - a^) ® 1® D ~ C 
and we represent the resolvents with the help of a new parameter 

D 



R{a 



(<) 



Jo 



da e~ a Q^)e 

c=l 



^■^=T^ c -^ iM ) 



(99) 



Note that, denoting Mf, q the minor of M with the line p and column q deleted, we have 



J-trlnM= ~J-lndetM 



d^r q - ~ m^ ~ ~ diFM4^ detM = diFM(- 1 ) P+9M P ? = \ M ~ l U ^nce 

the derivatives of the vertices are 



da 



trln[H(a«)l = ^trlnM(cr«) = - V ^f -^-tr InM(tr») 



ab 



da 



ab 






A 



1 ^(^a^ II <W') R ( a ®) 



N D 

np 

A ' 00 



p" 



jyD-1 



da e~ Q e VF 






,i_( (T M C _ cr (i)ct) 



ab 



, °v^^^ (i)c '- CT(i)c,t ) 



<9cr, 



(Oct 

bu 



trln[.R(a (i) )] 



iV^-i 



da e a 



-"^&^ )c -* (i)ct ) 



6a 
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JJ T r[ e - a \/^^ (CT(S>c '-' jWc ' t) 



(100) 



c'^c 



It follows that a tree made of one line of color c connecting two non ciliated vertices i 
and j will yield a contribution 



-A) 



daida 2 e - a i- Q2 Tr 



N ^ 1 / 

xnTr[e- Ql ^ lff( ' V - l ' ,c,,) ]nTr[e 



- ai \Z^ (ff<,)c ^ l!W) p~ a2 ^ (,7Wc - ff0W) 



-«2 A /^ir-r( CT 



h f^C3')j_„(j)jt 



c' — cr c 



(^^C 



(/^C 



, (101) 



which reproduces the equation (|98h for the tree with two vertices labelled i and j 
connected by one line. Note that this tree has 2D — 2 internal faces of colors d -^ c 
and one internal face of color c. On the other hand, a ciliated vertex writes as 



(J|fl((j( 



*)v)= r°da e -^n[ e_a ^^ (<r(il)c " CTCii)ct) ] c . ^w 

^° pn c=l P n 



reproducing eq. (f98|) for the tree with a unique ciliated vertex i\. 

The proof proceeds by induction on the number of vertices. Note that any tree can 
be obtained by adding one by one its edges. At each step two trees are joined by the 
new edges. As the lemma holds for the initial trees (as they have less vertices), when 
evaluating the derivative with respect to a^> and a^n one obtains a sum over terms, 
one for each occurrence of cr'** 1 and a^'\ i.e. a sum over all possible ways to join the 
two plane trees together into a plane tree with n vertices. Lemma [2] shows that the 
right hand side of equation (J98l) is reproduced. 

□ 



4.2.2 Proof of the second constructive expansion theorem [4] 

We add D fictitious integral over the unitary group U(N), i.e. we write 

W(J, J; X,N)= [ [dU l ] . . . f[dU D ] W(J, J; A, N) , 

Ju(N) J 



(103) 



which of course holds as J u(N ,[dU] = 1. Now, for all fixed U c , we perform the change 

of variables of Jacobian 1, a^ c —> U c ^a^' c U c . in eq. (|97p . The Gaussian measure 
invariant under this change of variables, hence 



is 



jf(n 



(i,j)eT n 



du lJ j / d/j, w , 
-> 

1 l r — 


> 

n n , / 


Z^ kl iV(fc+n-l)(D-l) 
k=0 i 

.. 2n-2+fc 

U )i®u(a) [ Yl da 

J q=l 


Z^ Z^ 

1 ,i 2 ...,i fc = l 7 - n O 



X 



/cgjrint^Oj 9 e<?(/ c ) 



(104) 
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k D 

X E Ti( J p\,-P? J n\,...n° Il U pfue( qi ) U h<lt cm )nl m 
pn,uv 1=1 c=l 



u c ( qi )v c (q ic(l) ) 



where the indices u c (qi) and v c (q^m) are summed. The integral over the unitary group 
of a product of matrix elements is, according to eq. (|34|) (see [Ml E5] for details) 






A 



E w sW TcCJ C ') II ^(^(o)- c (^(o)°-L(o^ c (o 
z=i 



0"c,T c 



E Wg(iV ' TcCJc ^ II ^ c fe)« c (9r c( o) 5 »: 



/=1 



i*Vc(0 



(105) 



where cr c and r c run over all the permutations of k elements (and in the second line we 
shifted both a c and r c by the permutation Cc 1 )- We obtain 



n n 

W(J, J; A. W) - £ -,(-*)-' E g yo^ixn-i) E E 

ra>l fc=0 'i«->'t=iT„ ,, 

2n-2+fc 



(106) 
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/( n dw<i ) [ d »wv( Tn ,u)i®°(<T) I ( n da e E '= T2+fe ° 9 
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x n t 4 n e- a <^^ q)c - i{qM) 
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x e (nn^tei^^ioi 

pn,uv 1=1 c=l 1=1 
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n 



Il{ J p},...pP J nl,...n?]l$ntpl cQ 
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TT e -«^(^ )c -"" )ct ) 



/C;i ( 9 ,)^ I ( 9Cc(0 ) 65oxt(7 -o^ geg(/ =i«(«!)^<(?e cC i))) 



« c fe)" c (95 c (o) 



The external sources group into trace invariants: each permutation a is one to one 
with a -D-colored graph Bg with Ik labelled vertices. The external strands recompose 
along the faces of the plane tree with external edges T^ L ~ to yield 



n n 

W(J,J;X,N) -E^-r'E ^^i),-!) E E 

71>1 ' fc = 

x E(n w g( Ar ' r ^" 1 )) ^m 



n! • ^ k\ AT(fc+n-l)(D-l) 

n>l fe=0 •i,i a ...,i fc =i7<3 i 



c=\ 
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2n-2+fc 



/ ( II du ) / d ^J{T n ,u)i®D{cr) I ( JJ da, 

• y ° (i.i)€T„ J J a=l 



e ^9 



E2n — 2+fc 
g = l Q 9 



It follows that the cumulants of the measure ^' are sums over graphs B 

Q(2k) 



(107) 



K \"-pu Inn • • • ^Pfc> ^n k ) 



OJp 1 uJn 1 . . . OJp k OJfl k 
C{B) 

B, fe(B)=fe P=l 



W( J, J; A, AT) 



j=j=o 



(108) 



where, denoting the (unique) permutation a associated to the L>-colored graph B with 
labelled vertices by <?(£>), we have 



1 i n 



(109) 



n>k 






!(f[Wg(iV,r c ^ 1 (fi))) /" ( [J <*«") j dn wii{ T n ,u)i^{a) 

-2 + fc -4 / ; 

da q ) e" E 9=! Q(J Tr e 9 V P=r l * 



2n-2+fc 

n 

9=1 

and Theorem [J] holds. 



FWfi*) qeq(fc) 



D 



4.2.3 Proof of the absolute convergence theorem [5] 

The operator a l ( q > c — a l y q > c ^ is an anti-Hermitian operator hence, denoting || • || 

I _ a /~^ ( ( ji(g)c_ (r «(g)ct) 

operator norm, e V Af 0-1 < 1 and 

Tr[ f[ e-^V^^( CTl(,)c - ff!(9)ct )j <iv|| J] e -^yf=§^{^"-^ M ) 



the 



</e?(/ c ) 



geq(f c ) 



< N 

q&q(f c ) 



n/^r( CT,(9)c - ffl(9)ct ) 



< JV 



(110) 



The integrals over a q and u lJ are bounded by 1 as well as the Gaussian integral (as 
(J> w i}(T n ,u)i® D ( (7 ) 1S normalized and positive). The Weingarten function is bounded, 
from eq. ([36]), by Wg(iV,r c o-~ 1 (i3)) < 2 — zrr— and we get a bound 

1 1 U 

w».^)i sE^i-' ^-V,) E EE 



n>k 



t 1 .t i ...,% h =lj^ l 
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xkl2 2Dk N ~2Dk+J2 c C(r c a- 1 (B))+\T(T^)\ ^^ 

The total scaling with TV is therefore, using eq. (12211 , 

-(k + n- 1)(D -l)-2Dk + Y, COcO-^B)) 

c 

+D + (n-l)(D-l)-J2 C(Cc) + J2 C{ - T ^c) 

c c 

= D - 2Dk - k(D - 1) + J2 C^a- 1 (B)) - £ C(&) + ^ C(r c £ c ) , (112) 

c c c 

which is bounded from Lemma [1] by 

D - 2Dk - k(D - 1) + (D + l)fc - C(B*) = D-2(D- l)k - C(B 3 ) . (113) 
Using the bounds in eq. (|73j) and eq. (j78j) we get 

1^(6,^)1 < ]V D -2p-l)fc-C(B ? ) fe! + l 2 2Dfc 

Y \\\ n - 1 (2D) n - \ (2 !!,| A: ~ 3)! „ • (114) 

Y V ; (n-fc)!(n + fc-l)! V ; 

A (not tight) bound on the combinatorial factor is 



(2n + k-3)l Wo _ , ,_ QNfe _ 2 (2n-i; 



. < (2n + k - 3) fc " V ;"•; ^ -- < (3?i) fc " 2 2 2n , (115) 
(n- fc)!(n + fc-l)! ~ v ; (n - fe)!(n + fc - 1)! ~ 

and we get 

|J!(B, /i$)| < AT^-^^- 1 )-C(B) fc p+i 2 2i3fe3fe-2 J- n^SDXr 1 . (116) 

n>fc 

The series is absolutely convergent for |A| < 2~ 3 L>~ 1 and the cumulants are bounded 
by 

\&(B,$)\ < N D-^D-1)-C(B)\ X \^K(B) , (117) 

for some constant K(B) independent of N (and independent of A for |A| small enough). 

□ 

4.2.4 Proof of corollary [1] 

Consider first the simpler problem of proving the function of a complex variable z = re tip 
defined as 

f(z) = i e- l 2dae- ae ~^- tarl/2x (118) 

Jo 

is analytic for ip E (— 7r,7r). First it is easy to see that / is bounded as \f(z)\ < — K?- < 

cos 2 

oo. Second, / respects the Cauchy-Riemann equations 

° i, \ f°° j -lie -ae"'7 1 ^ -jar 1 / 2 ! 

r Trf{z)= l dote 1 1 e -a—e 

or J 2 da 
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dae~^- — [ae~ c ** ^-wr 1 ' 11 * 



,£ 1 d 
2d, 

2 



dae- 4 t-(e- Qe Y 



, fe _ ae - ■) e - m v** = d f (119) 

ay 



hence it is analytic. For our case similar partial integration with respect to a shows 
that each function 

/2n-2+fe 

0=1 



^2n — 2+fc _— i 



x n Tr[ H e- Q '^/^^ (CT^(9,c - CTl(9)ct 



(120) 



respects the Cauchy-Riemann equations and, being bounded, is analytic in A. It follows 
that the full cumulant 



i i n 



121' 



n>k 






which is a sum of products of analytic functions, is analytic whenever the sum over n 
converges. Reproducing step by step the proof of Theorem [5] but taking into account 
that the integrals over a are bounded by — m- instead of 1, we obtain that the series 

° J cos -|- ' 

in eq. ([12ip is absolutely convergent for 



lAKCoos^V^" 1 



4.2.5 Proof of the mixed expansion theorem [6] 

Consider the contribution of a tree with external edges 



;i22) 
□ 



T E (T^J = N -D+2k(D-l)+C(B)-(k+n-D(D-l) fc , (j[ W g(iV, T c a~\B)) 

c=l 

X / ( II du J / d Vwv(T n ,u)l®n(a) 

/2n-2+k 
( e~ l zda q ) e~^"i= 



^2n — 2+fc _— t 



e "2"o„ 



x J] Tr[ f[ e -^V^^^ (9)c -^ M ) 
We use a Taylor expansion 



(123) 



1 Yd q 



/(^e^) = £ [_ /(v /^e 



g=0 



t=0 
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1 



d s 



+ (I^/o (1 - tr V 



f(^/m^))dt. 



(124) 



We must evaluate the derivative with respect to t acting on T (T CJ -). In each operator 

y|A| multiplies a difference a — o^ and the derivative with respect to t acting on an 
exponential computes to 



d 



n Tr f n e^ v V^& i(9)c - CTi(g)ct )" 

e n Tr [ n e^^^ (ff,(, ^ ff,(,M: 

1 V-^/ i(q) c 8 



"g£g(/ c ) 



^ 2t 2-*> \°^ 

qeq{n ab 



i(q)c + a ba 



ab 



da, 



*(9)ct 

ba 



;Tr[ TT e^^^ (ff,( "" ff '" M) 



ge 9 (/ c ) 



n D 



EE^E( ct 

i=l c=l ab 



r (0= 



s 



Oct 



+ ^ 



(0=t 



s 



'ab „ (i) e ^"ba (<) ct 



aft 



<9a 



/)« 



n T i n e 

Integrating by parts the Gaussian integral we get 
d 



^vVi^M ffl(9)c - CTl(9)ct ) 



(125) 



dfi % 



It J ' 

dV W i3(T n ,u)l® D \ 



u )i*u(a) || Tr^ || e V n"-^ 
/ c e^"(r n ° l? ) <?6g(/ c ) 



n D 



X 



c=l a fe v 0"7 ba C^ ab 



d d d d 

+ 



li)c\ fsjfic 



d^ c da o r 



n Tr[ n e -«^/A(" i(,)c - ff,( ' M ) 

f c ^(Tf, t) q(iq{ -^ 



(126) 



The action of the derivatives on a product of traces has been evaluated in lemma [U 
The sums over i, j and c yields a sum over all the possible ways to add a colored loop 
edge to the tree TP *. The two new auxiliary parameters f3± and 71 (which we relabel 

as two supplementary a parameters) don't have any e _J 2" in the measure, hence we 
need to explicitly add them. Thus we get 



1 , t\\\e t ^\ I 



JUi.c 



2n-2+k+2 

n 



_ xp2n-2+fc+2 -t 

e *da„ e ^^i 
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II H II e"^^^ (ff!(9)c "^ )ct) ] . (127) 

Note that both derivative terms lead to the same contribution (the loop edge does not 
have any orientation) hence 



d_ 

7ft 



D 



T E (T^J = N -D+2k(D-l)+C(B)-(k+n-l)(D-l) fcl JJ W g(iV, T^ 1 ^) 



X 



/ ( II du ) I d ^wH{T n ,u)i^{a) yi 
(i,j)&r n J c={(r AI 

„ 2n-2+fc+2 

/( ]1 e^da q )e-^ 2+k+2 ^° 



c=l 



£={(hJA)Y,ci 



N D- 



T )w^(T n ,u) 



X J] Tr[ f[ e -"*Styf^(o*M°-e«M) 

and the derivative of order s is 



(128) 



lt : 



D 



T E (T^J = N -D+2k(D-l)+C(B)-(k+n-l)(D-l) jfcl ( JJ W g(iV, TcCr" 1 ^) 



X 



/ ( n du j / d /*tt«(T B) u)i® D ( ff ) x] 

2+A 

n 



c=l 



(*J')6T, 
3 _ 2n-2+fc+2s 

! / 



v On^^c^, 



£;ci...c s p=l 



M 



9=1 



da 9 ) e ^9=1 



x n Tr[ n e-w^**-*^ 



(129) 



where £ = { (ji , j[ ) , . . . (is, is)} runs over all the possible ways to decorate T^ L ? with 
unoriented loop edges C and a, . . . c s run over the possible colorings of the 2s loop 
edges. Taking into account that the Gaussian measures are normalized we evaluate 



T E (T° 
Idtt K "' t >'"Jt= 



D 



N ^ D+2k(D ^ )+c{ B)-(k + n-i)(D-i) kl rY[Wg(N,T c a-\B) 



c=l 

o 






(i,j)&T n 



C;ci...c q p=l 



(-^)Mn wg ( iv ' rcac ~ 1(B) )) e A n d« y )n«^»(r B1 u) 

c=l C,ci...c n J ° (i.i)GT„ p=l 



Ar -D+2fc(D-l)+C(S)-(fe+r l -l)(D-l)+|^(r n » > - £ )| ^ /^ 

which yields the terms r«'(TP ~) and T EL >( q > (7~ -, ,.) in Theorem EJ The rest term is 






dt(l - t) 



s-1 



^y- Z)+2Jfc(Z)-l)+C(B)-(fe-|-n-l)(Z)-l) 



38 



D 



k\(llWg(N,T c a-\B))) Yl 



c=\ 



C;ci...c s 



A 



N D-1 



( II dv? j }f[w>*>p(T n ,u) fdfi 

(i,j)£T n P =l 



w i i{T n ,u)ie> D ( (7 ) 



oo / 2n-2+k+2s 
9=1 



2n-2+fc+2s -i* 
9 = 1 e * a 1 






[131) 



which establishes the mixed expansion. Concerning the bounds, we first bound the 
termT(?)(7^) 



r (9) (0= £ ^ 



£,Cl...Cq 



£> 



N D 



c=l 



(132) 



(M)eT„ p=i 

We bound the product of Weingarten functions by 2 2Dk N~ 2Dk+c( - Tc(Tc ^\ The inte- 
grals over u 13 are bounded by 1, hence 



\ T (q) (T l _ ^ J.ln2Dk V^ 

I 1 \'n,<„f)\ r ,\ AT(D~l)q K - Z 2^ 



(133) 



£,ci...c„ 



gl N(D-l)q 

xN -D+2k{D-l)+C(B)-2Dk+C(T c a- 1 {B))-{k+n-l){D-l) + \HTZ,r,c)\ . 

By eq. (|3ip we find that the scaling with N is bounded by 

-D + 2k(D - 1) + C{B) - 2Dk + C{T c a~ l {B)) - (k + n - 1)(D - 1) (134) 

+D + (n-l)(D-l)-J2 C(&) + Y, C ( T ^ + 9 - «P - 1) < -?P " 2 ) . 

c c 

where we used lemma [TJ Thus 

1X19 „„, (Or, 4-9«4- fc — 3"\! 

(135) 



^n,*,WI-= ArafD-2^" 4 ^ g!(2n + £;-3)! ' 



Ni( D ~ 2 ) 

For the rest term, R( S '(T CJ -) we use similar bounds and taking into account that the 
integrals over a are bounded by — % and the integral over t is bounded by 1 we find 



\R {S) (TZ*)\ < 



i 



^\2n+2s+k-2 jys(D-2) 



(cos f ) 



|A| S fc!2 2^^ (2n + 2 S + fc-3)! 



(s-l)!(2n + fc-3)! 



(136) 
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4.2.6 Proof of the Borel summability theorem \7\ 

We now show that the rescaled cumulants 

K(B,N) = 7V- D+2fc ( D - 1 )+ c '( B )^(B ) 4 ) ) 

= E^(-^ 1 E EE r *(0 



n>k 



111 



(137) 



i 1 ,i 2 ...,i k =lj£> i T 



are Borel summable in A uniformly in N. First, the corollary Q] ensures that the series 
(f443"|) is absolutely convergent for |A| < (cos ^) 2~ 3 D~ l , hence it certainly is absolutely 
convergent in a Borel disk of fixed radius. 

Second, using the mixed expansion theorem we write 



jw^E^-a)"- 1 E EE[E T(9) (0+* (s) (^ 



n>k 



»l.»2— .»fc =1 7^ t ? "? =1 



(138) 



We perform a Taylor expansion of K(B, N) in A up to order r > k. All the terms 
corresponding to trees with n > r + 1 are in the reminder, hence for them we use the 
mixed expansion for s = 0. For the terms corresponding to trees with n < r + 1 we use 
the mixed expansion up to order s = r — (n — 1). The explicit terms T^ q '(l~° S) yield 
a series in A which is nothing but the Taylor expansion of the rescaled cumulant. The 
reminder term of the Taylor expansion of K(B, N) writes as 



r .. n 

w,a) =E^(- A ) n - 1 E EE« (r+1 - n) (^ 



-) 



n=k 



i K i 2-< i k= 1 Tn,i. T 
id^id' 
n 



+ E ^r- 1 E ££* (0) cC 



(139) 



n>r+l 



i lt i 2 . .i k =lj2 b T 
l d¥=i d ' 



The terms with n > r + 1 then admit a bound 

I E h^r- 1 £ EE* (0) C7£, 



n>r+l 



»! 






<^ ( n _ fc)!( n + fc _ i)! ( CQS <p^n-2+k K - 

1 



2Dfc 



n>r+l 

< k \ D +^™ y: \M n -\m n -\ 

n>r+l [ cos 2) 

|A| ^ 



(3n) 



fc-2 2ra 



(cos ^) ^ (cos ^) 



(140) 



for some constant K and |A| < 2 3 D 1 (cos i0 . The terms with r < n admit a bound 

r - n 

lE^r 1 E E^-'cfo,: 



n=k 
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r 



< 



^|Ar" 1 (2D) 



n _ 1 (2n + k-3)\ 



(n - k)l(n + k - 1)1 

n=k 

Uir+1-njy- (r+l-n)(D-2)^,|C+l 2 2Dfc 



(cos |) 



2n-2+k+2(r+l-n) 



r+ i- w [2n + A; + 2(r + l-n)-3]! 

(r-n)![2n + lfe-3]! l J 



Taking into account that 

(2n + fe-3)! 



< S 2 "^- 1 ^! < 3 3r fc! 



(n- k)\(n + k- 1)\ 

< 3 2r+fc-i (r _ n + 2)! < 3 3 > + 1)! (142) 



(r - n)l[2n + k - 3]\ 
these terms are bounded by 

|A| r 



(cos f ) 

Overall we thus derive a bound 

L\ |r 
(cos f ) 

for some constants K and <t which proves the theorem [7J 



__ k[ D+2 2 2Dk D r 3 6r^ + jjjgr ^43) 



i?jv,,(B, A) < K - \; s2r+k (T r r\ , (144) 



□ 



4.2.7 Proof of the 1/N expansion theorem [8] 

Using the mixed expansion up to order s for every T^ L * leads to the rest term 

1 n 

R${BA) = Y.h { - x)n ~ 1 E EE* (s) (0' ( 145 ) 

n>fc i 1 ,i 2 ...,i fe =l7<i t r 

hence using the bound on RS S > (T^ ^) we get 

n 

i^c^A^E^r 1 E EE 

n>fc i 1 ,i 2 ...,i k ^l-f-0 L f 

x 1 |A|- m 2Dk DS) (2n + 2, + fc-3)! 

(cos |) 2 "+ 2s + fc " 2 iV-P-2) ^ J (s - l)!(2n + k - 3)! ' l j 

By the combinatorial countings eq. (f73j) we get 



|flW(B,A)| < -^ T (A;! D + 1 2 2Bfe D s )^— , Vf^LV +n ' 

" JI "iV^ 2 ) 1 (cosf)^(cosfr 

,, f2 ™-i (2n + fc-3)! (2n + 2, + fc-3)! 

V ' (n-k)\(n + k-l)\(s-l)\(2n + k-3)\ ,{ ' 
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that is, denoting K and a two constants depending only on k, 



N s(D-2) " / £ x2.+3k-l 



COS ^ ] 



v / |2DA| y (2q + 2s + 3k-3)\ 



^cosf) 2 ^ (s-l)\q\(q + 2k-l 



A non optimal bound on the combinatorial factor is 

(2g + 2s + 3/g-3)! 2 g+2s+3fc-3 / 14g x 

( S -l)\q\(q + 2k-l)\- {S + L) - k - b ' li4yj 

and, as the sum over q converges absolutely for |A| < 5 _2 2~ 1 L>~ 1 (cos ^) , we get 

|i$(S,A)| < ^^' ^ ( L 2 ) (co ^ + 2 ^-i • ( 15 °) 

for some constants K and cr. 

□ 
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